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Abstract. We present a complete recipe to extract the density-density correlations 
and the static structure factor of a two-dimensional (2D) atomic quantum gas from 
in situ imaging. Using images of non-interacting thermal gases, we characterize and 
remove the systematic contributions of imaging aberrations to the measured density- 
density correlations of atomic samples. We determine the static structure factor and 
report results on weakly interacting 2D Bose g ases, as well as strongly interacting 
gases in a 2D optical lattice. In the strongly interacting regime, we observe a strong 
suppression of the static structure factor at long wavelengths. 
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1. Introduction 

Fluctuations and correlations result from the transient dynamics of a many-body system 
deviating away from its equilibrium state. Generally, fluctuations are stronger at 
higher temperatures and when the system is more susceptible to the external forces 
(as governed by the fluctuation-dissipation theorem, see ^E])- Local fluctuations and 
their correlations can thus be a powerful tool to probe thermodynamics, and to identify 
phase transition of a many-body system due to the sudden change of the susceptibility 
to the thermodynamic forces. 

Measurement of fluctuations and correlations on degenerate atomic gases can reveal 
much information about their quantum nature [3]. Experiment examples include the 
quantum statistics of the atoms [H El O [3 |8], pairing correlations [9j and quantum 
phases in reduced dimensions [1^ In these experiments, images of the sample are 
taken after the time-of-flight expansion in free space, from which the momentum-space 
correlations are extracted. 

In situ imaging provides a new and powerful tool to examine the density fluctuations 
in real space [121 [El HH [IS [IS [El [El [El [2^, offering a complimentary description 
of the quantum state. This new tool has been used to resolve spatially separated 
thermodynamic phases in inhomogeneous samples. From in situ measurements, both 
Mott insulator density plateaus and a reduction of local density fluctuations were 
observed [T3l [El [El [21]. Furthermore, a universal scaling behavior was observed in 
the density fluctuations of 2D Bose gases [22] . 

Precise measurements of spatial correlations, however, present signiflcant technical 
challenges. In in situ imaging, one typically divides the density images into small unit 
cells or pixels and then evaluates the statistical correlation of the signals in the cells. 
If both the dimension of the cell and the imaging resolution are much smaller than the 
correlation length of the sample, the interpretation of the result is straightforward. In 
practice, because the correlation length of quantum gases is typically on the order of 
1 //m, comparable to the optical wavelength that limits the image resolution, interpreting 
experimental data is often more diflicult. Finite image resolution, due to either 
diffraction, aberrations or both, contributes to systematic errors and uncertainties in 
the fluctuation and correlation measurements. 

In this paper, we present a general method to determine density-density correlations 
and static structure factors of quantum gases by carefully investigating and removing 
systematics due to imaging imperfections. In Section 2, we review the static structure 
factor and its relation to the real space density fluctuations. In Section 3, we describe 
how the density fluctuation power spectrum of a non-interacting thermal gas can be used 
to calibrate systematics in an imperfect imaging system, and show that the measurement 
can be explained by aberration theory. In Section 4, we present measurements of density 
fluctuations in weakly interacting 2D Bose gases and strongly interacting gases in a 
2D optical lattice, and extract their static structure factors from the density-density 
correlations. 
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2. The density-density correlation function and the static structure factor 

We start by considering a 2D, homogeneous sample at a mean density n. The density- 
density correlation depends on the separation ri — r2 between two points, and the static 
correlation function iy{r) is defined as p3] 

z/(ri - r2) = n~^{5n{ri)5n{r2)) 

= 5(ri - rs) + n-^(^^t(r,)^t(r2)^r(ri)^r(r2)) - n (1) 

where (...) denotes the ensemble average, and 5n{r) = n(r) — n is the local density 
fiuctuation around its mean value n. The Dirac delta function 5{ri — r2) represents the 
autocorrelation of individual atoms, and (^'''(ri)^^'''(r2)^^(ri)^(r2)) = G^^^(ri — r2) is 
the second-order correlation function [24]. When the sample is completely uncorrelated, 
only atomic shot noise is present and z/(ri — r2) = 5(ri — r2). At sufficiently high phase 
space density, when the inter-particle separation becomes comparable to the thermal de 
Broglie wavelength A^^ or the healing length, density-density correlation becomes non- 
zero near this characteristic correlation length scale and iy{r) deviates from the simple 
shot noise behavior. 

The static structure factor is the Fourier transform of the static correlation function 
[231 ES] 

S{k) = J zy(r)e-^*^"rfr, (2) 

where k is the spatial frequency wave vector. We can rewrite the static structure factor 
in terms of the density ffuctuation in the reciprocal space as 



_ {5n{k)5n{-k)) _ {\5n{k)\') 

- ^ - , (<^J 

where 5n(k.) = J 5n{r)e~^^'^dr^ and is the total particle number. Here, 5n(— k) = 
5n*(k) since the density ffuctuation Sn{r) is real. The static structure factor is therefore 
equal to the density ffuctuation power spectrum, normalized to the total particle number 
N. A non-correlated gas possesses a structureless, ffat spectrum S(k) = 1 while 
a correlated gas shows a non-trivial S'(k) for k near or smaller than inverse of the 
correlation length 

The static structure factor reveals essential information on the collective and the 
statistical behavior of thermodynamic phases [2[ [25l [26l [27] . It has been shown that, 
through the generalized fluctuation-dissipation theorem [28], the static structure factor 
of a Bose condensate is directly related to the elementary excitation energy e{k.) as 
[251 [29] 

where m is the atomic mass, T is the temperature, h is the Planck constant h divided by 
27r, and is the Boltzmann constant. See references [25l[26l[27] on the static structure 
factor of a general system with complex dynamic density response in the frequency 
domain. 
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Figure 1. A comparison between physical length scales and measurement length 
scales. In (a), each atom is represented by a color circle with its diameter equal to the 
thermal de Broglie wavelength, A^^^. ^ is the correlation length. An ideal measurement 
detects atom with perfect resolution, (b) shows the experimental condition where the 
image resolution is larger than the other length scales. Here, the image of an atom 
forms a resolution limited spot (dashed circles), and is large compared to the correlated 
area and the CCD pixel area. The grid lines represent the CCD pixel array. 



Previous experimental determinations of the static structure factor in the zero- 
temperature limit, based on two-photon Bragg spectroscopy, have been reported for 
weakly interacting Bose gases j30l EI] and strongly interacting Fermi gases [32]. Here, 
we show that ^^(k) at finite temperatures can be directly determined from in situ density 
fluctuation and correlation measurements. 

Experimental determination of ^^(k) from density fluctuations is complicated by 
artiflcial length scales introduced by the measurement, including, for example, flnite 
image resolution and size of the pixels in the charge-coupled device (CCD). Figure 1 
shows a comparison between the measurement length scales (the resolution limited spot 
size and CCD pixel size), the correlation length ^, and the thermal de Broglie wavelength 
XdB- Ideally, a density measurement should count the atom number inside a detection 
cell (pixel) with sufliciently high image resolution, and the dimension of the cell should 
be small compared to the atomic correlation length. In our experiment, the image 
resolution is determined by the imaging beam wavelength A = 852 nm, the numerical 
aperture N.A. = 0.28, and the aberrations of the imaging system. The image of a single 
atom on the CCD chip would form an Airy-disk like pattern with a radius comparable 
to or larger than A^^^ or ^. The imaging magniflcation was chosen such that the CCD 
pixel size ^/A = 0.66 /xm in the object plane is small compared to the diffraction limited 
spot radius ^1.8 /xm. The atom number Nj recorded on the j-th CCD pixel is related 
to the atom number density / drn{r)V{Tj — r), assuming the point spread function V{t) 
is approximately flat over the length scale of a single pixel, 

riexpirj) = J dvn{v)V{v^ - r), (5) 

where Yj is the center position of the j-th pixel in the object plane, n(r) is the atom 
number density distribution, and the integration runs over the entire x — y coordinate 
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space. The atom number fluctuation measured at the j-th pixel is related to the density- 
density correlation as 

{5N^) ^A^JdrJ dv'{5n{v)5n{Y'))V{Y, - r)P(r, - r^, (6) 

where 5Nj = Nj — (Nj) is the atom number fluctuation around its mean value {Nj). 
In the Fourier space, Eq. (5) can be written as 

5ne.p(k0^5n(kOP(kO, (7) 

where SuexpO^i) = Zlj SNje~'^^^'^^ is the discrete Fourier transform of 5Nj^ approximating 
the continuous Fourier transform. Here, k/ = ^(Z^, ly)^ L is the linear size of the image, 
Ix and ly are integer indices in k-space. From Eq. (3) and (7), the power spectrum of 
the density fluctuation is related to the static structure factor as 

{\5ne.,{ki)\')^NS{ki)M\ki), (8) 

where the modulation transfer function A^(k) = |P(k)| accounts for the imaging 
system's sensitivity at a given spatial frequency k, and is determined by the point 
spread function. Also, from Eq. (6), the pixel- wise atom number fluctuation is related 
to the weighted static structure factor integrated over the k-space, 

{SNf) ^ I dkS{k)M'{k). (9) 

Generalization of the above calculations to arbitrary image resolution and detection 
cell size is straightforward. In addition to convolving with the point spread function, 
the measured atom number density must also be convolved with the detection cell 
geometry. Equation (5) can therefore be written as Nj/A = J rfkn(k)P(k)'H(k)e^*^'^^ , 
where 'H(k) = fj^e^^'^dr /A and the integration goes over the area A of the detection 
cell. This suggests simply replacing A^^(k) by A^^(k)'H^(k) to generalize Eq. (8) and 
(9). Finally, in all cases, the discrete Fourier transform deflned in Eq. (7) should well 
approximate the continuous Fourier transform for spatial frequencies smaller than the 
sampling frequency \j\J~A. 

We view the factor A4^(k)'H^(k) as the general imaging response function, 
describing how the imaging apparatus responds to density fluctuations occurring at 
various spatial frequencies. To extract the static structure factor from in situ density 
correlation measurements, one therefore needs to characterize the imaging response 
function at all spatial frequencies to high precision. Since our pixel-size is much smaller 
than the diffraction and aberration limited spot size, we will from here forward assume 
'H^(k) = 1; 'H^(k) decays around k ^ 4:/^/A = 6 fim~^ (1/e radius), which is much 
larger than k = 27rN.A./A = 2.1 //m~-^, where Ai^(k) terminates. 

3. Measuring the imaging response function A^^(k) 

In this section, we show how to use density fluctuations of thermal atomic gases to 
determine the imaging response function A^^(k). Other approaches based on imaging 
individual atoms can be found, for example, in Refs. [I33[ l3^ . 
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3.1. Experiment 



Measuring density fluctuations in low density thermal gases provides an easy way to 
precisely determine the imaging response function. An ideal thermal gas at low phase- 
space density has an almost constant static structure factor up to /c = A^^ [24] which, in 
our case, is larger than the sampling frequency \j\J~A. Therefore, the density fluctuation 
power spectrum derived from an ideal thermal gas reveals the square of the modulation 
transfer function, as indicated by Eq. (8). 

We prepare a 2D thermal gas by flrst loading a three-dimensional ^^^Cs Bose- 
Einstein condensate with 2 x 10"^ atoms into a 2D pancake-like optical potential with trap 
frequencies cj^ = 27r x 2000 Hz (vertical) and cj^ = 27r x 10 Hz (horizontal) [211 [22]. We 
then heat the sample by applying magnetic fleld pulses near a Feshbach resonance. After 
suflBcient thermalization time, we ramp the magnetic fleld to 17 G where the scattering 
length is nearly zero. The resulting thermal gas is non-interacting at a temperature 
T = 90 nK and its density distribution is then recorded through in situ absorption 
imaging [22j. 

We evaluate the density fluctuation power spectrum A^g^^(k/) = (|5nea:p(k/)p) using 
60 thermal gas images (size: 256 x 256 pixels). Figure 2(a) shows a sample of the noise 
recorded in the images. Outside the cloud (whose boundary roughly follows the dashed 
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Figure 2. Determination of the imaging response function A^^(k) from in situ images 
of 2D thermal gases, (a) Image noise of a 2D thermal gas. The dashed ellipse encircles 
the location of thermal atoms, (b) Noise power spectrum evaluated from 60 images, 
using discrete Fourier transform defined following Eq. (7). Zero spatial frequency is 
shifted to the image center, (c) Fit to the image noise power spectrum using imaging 



response function defined in Eq. (12) and aberration parameters defined in Eqs. (10) 



and ( pTj ). (d) Sample line-cuts of experiment (circles) and fit (solid lines), with cutting 
angle indicated in the graph. The profiles are plotted with offset for clarity. Image 
size: = 256 x 256 pixels. 
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line), the noise is dominated by the optical shot noise, and is therefore uncorrelated 
and independent of spatial frequency. In the presence of thermal atoms, we observe 
excess noise due to fluctuations in the thermal atom density. The noise power spectrum 
is shown as an image in Fig. 2(b), with line-cuts shown in Fig. 2(d). We note that 
the power spectrum acquires a flat offset extending to the highest spatial frequency, 
due to the photon shot noise in the imaging beam. Above the offset, the contribution 
from atomic density fluctuations is non-uniform and has a hard edge corresponding to 
the flnite range of the imaging response function. Close to the edge, ripples in the 
noise power spectrum appears because of aberrations of the imaging optics, discussed 
in later paragraphs. Finally, the bright peak at the center corresponds to the large scale 
density variation due to the flnite extent of the trapped atoms, and is masked out in 
our following analysis. 

To fully understand the imaging response function with imaging imperfections, we 
compare our result with calculations based on Fraunhofer diffraction and aberration 
theory [35j as described in the following paragraphs. 

3.2. Point spread function in absorption imaging 

We consider a single atom illuminated by an imaging beam, the latter is assumed to be 
a plane wave with a constant phase across both the object and the image planes. The 
atom, driven by the imaging electric fleld, scatters a spherical wave (dark fleld) which 
interferes destructively with the incident plane wave [36j. The dark fleld is clipped by 
the limiting aperture of the imaging optics and is distorted by the imaging aberrations. 
An exit pupil function p is used to describe the aberrated dark fleld at the exit of the 
imaging optics [35j, and its Fourier transform p(k.) with k oc R describes the dark fleld 
distribution on the CCD chip, where R is the position in the image plane. The image 
of a single atom is then an absorptive feature formed by the interference between the 
dark fleld and the incident plane wave in the image plane. 

We extend this to absorption imaging of many atoms with density n(r) = ^iS{r — 
r^), where is the location of the i-th atom in the object plane. The total scattered 
fleld in the image plane i^AE" = 6]9(k — k^), with each atom contributing a dark fleld 
amplitude e, and k relates to the position r in the object plane through k = r/arf, where 
a is the radius of the limiting aperture and d = A/(27rN.A.). The dark fleld e oc e^^^E^ 
is proportional to the incident fleld £"0, and carries with a phase shift associated with 
the laser beam detuning from atomic resonance. For a thin sample illuminated by an 
incident beam with intensity /q, the beam transmission is = |£'o + A£'p/|£'op ^ 
1 + 2^[/\E/E^] and the atomic density n^xp oc - ln(t2) + (1 - t^)Io/Isat [SZl |22] leads to 
riexp oc -2{l + Io/IsaM^E/Eo] oc Ei^e'^^p{k-ki)]. Here, refers to the real part 
and I sat is the saturation intensity for the imaging transition. Comparing the above 
expression to Eq. ([5]), we derive the point spread function as V{r) oc ^[e'^^^p(k)]\\^=r/adi 

i We consider the density n of the 2D gas in a range that the photon scattering cross section remains 
density-independent [38^ . 
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in contrast to the form |p(k)p in the case of fluorescence or incoherent imaging. 

When the dark fleld passes through aberrated optics, neither the amphtude nor 
the phase at the exit pupil is uniform, but is distorted by imperfections of the imaging 
system. To account for attenuation and phase distortion. We can modify the exit pupil 
function as 

p(r„^,) = [/(r,/a,^,)e^®(^-/«'^-), (10) 

where and 9p are polar coordinates on the exit pupil, U (p, 9) is the transmittance 
function, and Q{p^9) is the wavefront aberration function. We assume U to be 
azimuthally symmetric and model it as U{p) = H{1 — p)e~^^/'^^, where H{x) is the 
Heaviside step function setting a sharp cutoff when reaches the radius a of the limiting 
aperture. The factor e"^^/'^^, with 1/e radius = ar, is used to model the weaker 
transmittance at large incident angle due to, e.g., flnite acceptance angle of optical 
coatings. For the commercial objective used in the experiment, we need only to include 
a few terms in the wavefront aberration function 

e(p, 9) ^ Sop^ + V cos{29 - 20) + (11) 

where the parameters used to quantify the aberrations are: Sq for spherical aberration, 
a for astigmatism (with (j) the azimuthal angle of the misaligned optical axis) , and /5 for 
defocusing due to atoms not in or leaving the focal plane during the imaging. 



Using the exit pupil function in Eq. (10), we can evaluate the point spread function 
via V{r) oc ^[e'^^^p{^)]\\i=r/ad with proper normalization. We can also calculate the 
modulation transfer function (k) = |P(k)| . In fact, determination of any one of ]9(r^), 
P(r), or A^(k) leads to a complete characterization of the imaging system including its 
imperfections. 

3.3. Modeling the imaging response function 



We flt the exit pupil function p in the form of Eq. (10), using a discrete Fourier transform, 
by comparing 

M%, = \TTme^''mp)W (12) 
to the thermal gas noise power spectrum Ail^p shown in Fig. 2(b). Here, J-'T{.) denotes 
the discrete Fourier transform. Figure 2(c) shows the best flt to the measurement, 
which captures most of the relevant features in the experiment data. Sample line-cuts 
with uniform angular spacing are shown in Fig. 2(d). This experimental method can in 
principle be applied to general coherent imaging systems, provided the signal-to-noise 
ratio of the power spectrum image is sufliciently good to resolve all features contributed 
by the aberrations. Moreover, one can obtain analytic expressions for the point spread 
function and the modulation transfer function once the exit pupil function is known (see 
Appendix A). 

Having determined the imaging response function, one can remove systematic 
contributions from imaging imperfections to the static structure factor as extracted 
from the power spectrum of atomic density fluctuations, see Eq.(8). 
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Figure 3. Illustration of the patch selected for the static structure factor analysis, 
(a) shows a typical cloud image of 200 x 200 pixels. The selected patch is located at 
the center of the cloud, bounded by a box with an area of 32 x 32 pixels, (b) shows 
the density fluctuations inside the patch. 

4. Measuring density-density correlations and static structure factors of 
interacting 2D Bose gases 

We measure the density-density correlations of interacting 2D Bose gases based on 
the method presented in the previous sections. This study is partially motivated by 
a finding in our earlier work that the local density fluctuation of a 2D Bose gas is 
suppressed when it enters the Berezinskii-Kosterhtz-Thouless (BKT) fluctuation and the 
superfluid regions [22j . We attributed this phenomenon to the emergence of long density- 
density correlation length exceeding the size of the imaging cell and the resolution. This 
results in a smaller pixel- wise fluctuation 6N'^/A than the simple product of the thermal 
energy /c^T and the compressibility as is expected from the classical fluctuation- 
dissipation theorem (FDT) [Ij. Below, we present a careful analysis of the density- 
density correlations of interacting 2D Bose gases and discuss the role of correlations in 
the FDT. 

To extract local properties from a trapped sample, we limit our analysis to a small 
central area of the sample where the density is nearly flat. In addition, the area is 
chosen to be large enough to offer sufficient resolution in the Fourier space. We choose 
the patch size to be 32x32 pixels. Figure |3] shows a typical image and the density 
ffuctuations inside the patch. 

To ensure that we obtain an accurate static structure factor using the small patch, 
we perform a measurement on a non-interacting 2D thermal gas at a phase space 
density nX^^ = 0.5 and compare the measured static structure factor to the theory 
prediction [24j. We flrst calculate the imaging response function A^^(k) for a patch 
size of 32 x 32 pixels and divide the thermal gas noise power spectrum by Ai^(k). The 
resulting spectrum should represent the static structure factor of an ideal 2D thermal 
gas. In Fig. [4j we plot the azimuthally averaged static structure factor with data points 
uniformly spaced in fc, up to the resolution limited spatial frequency k = 27rN.A./A. The 
measured static structure factor is flat and agrees with the expected value of S{k) ^1.3 
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for k < X^^ = 2 /xm"^ 

Applying the same analysis to interacting 2D Bose gases, we observe very different 
strengths and length scales for the density fluctuations. In Fig. |4](a-c), we present the 
single-shot image noise of samples prepared under three different conditions: weakly 
interacting gases at the temperature T = 40 nK (below the BKT critical point), with 
dimensionless interaction strengthjf] 5^ = 0.05 and 0.26 (phase space density nX^^ = 9 
and 7); and a strongly interacting 2D gas at the temperature T = 8 nK, prepared in a 
2D optical lattice at a mean site occupancy number of 2.6, and a depth of 7 Er^ where 
Er = h X 1.3 kHz is the recoil energy. Due to the tight conflnement, the sample in 
the optical lattice has a high effective interaction strengthf g^^ = 1.0 [39j. Details on 
the preparation of the 2D gases in the bulk and in an optical lattice can be found in 
Refs. [22J and [39j, respectively. 

The difference in the density fluctuations shown in Fig. [4](a-c) can be characterized 
in their static structure factors shown in Fig. [4](d). We observe positive correlations 
above the shot noise level S{k) = 1 in the two weakly interacting samples. The one 
at g = 0.05 shows stronger density correlations at small k than does the sample at 
g = 0.26. The enhanced density correlations S{k) > 1 at low momenta are expected 
since the thermally induced phonon excitations can populate states with length scale 
1/k longer than the healing length ^ = If-^/ng. For gases with stronger interactions, 
excitations cost more energy and the excited states are less populated. At smaller g^ 
the correlation length is longer and, therefore, the static structure factor decays at a 
smaller k. 

The most intriguing observation is the negative correlations S{k) < 1 in the 
strongly interacting gas with g^^ = 1.0. We observe a below-shot-noise spectrum at low 
momentum fc, showing that long wave-length excitations are strongly suppressed due to 
a stronger interaction energy nh^gef^/m = /c^ x 34 nK compared to the thermal energy 
ksxS nK. As the momentum k increases, the excitation populations gradually return to 
the shot noise level. Our observation is consistent with the prediction in Ref. [29j that 
when the thermal energy drops below the interaction energy, global density fluctuations 
in a superfluid are suppressed. 

Finally, we discuss the contribution of flnite density-density correlations in the 
FDT. Including correlations, we can write the FDT as /c^T/^(r) = /(5n(r)5n(r'))rfr' = 
n(r)5'(0) [29j. We compare the measured static structure factor, extrapolated to 
zero-/c, to the value of fc^T/^/n, where = dn/dfi is the experimentally determined 
compressibility [22], and indeed flnd that nS'(O) equals to fc^T/^ to within our 

§ Following the calculation in Ref. [24], we find the static correlation function of an ideal 2D thermal 
gas is v{r) = 5{r) + \gi{z^e~^'^ 1)A^^, where z = e^^^^^ is the local fugacity and 

g^{x,y) = YlV=i^^y^^^ /^^ generalized Bose function. Fourier transforming v{r) to obtain 

the static structure factor S{k)^ we find S{k) ^12 remains flat for kX^B < 1- 

t The dimensionless interaction strength of a weakly interacting 2D Bose gas is ^ = ^/87^as/lz^ where 
as is the atomic scattering length and 1^ = \/h/muJz is the vertical harmonic oscillator length. For a 
2D gas in a 2D optical lattice, the effective interaction strength is ^eff = rriUP/h^ [39 , where U is the 
on-site interaction and / is the lattice constant. 
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Figure 4. Density fluctuations and the static structure factors of 2D Bose gases, (a) 
and (b) Image noise of weakly interacting 2D Bose gases in the superfluid phase at 
dimensionless interaction strength g = 0.05 and 0.26. (c) Image noise of a strongly 
interacting 2D Bose gas at ^eff = 1.0 prepared in a 2D optical lattice at a depth of 
7 Er. (d) shows the static structure factors extracted from the noise power spectra of 
interacting 2D gases as shown in (a) (black circles), (b) (red squares), and (c) (blue 
triangles). The static structure factor of an ideal thermal gas at phase space density 
nX^^ = 0.5 (open circles) and the expected value of S{k) ^1.3 (gray dashed line) 
are plotted for comparison. Solid lines are the guides to the eye. Vertical dashed line 
indicates the resolution limited spatial frequency k = 27rN.A./A = 2.1 /im~^. 



experimental uncertainties of 10 ^ 20% for all three interacting samples. This agreement 
shows that the measured correlations and thus the static structure factor can be linked 
to the thermodynamic quantities via the FDT. Our ability to determine S{0) and k 
from in situ images also suggests a new scheme to determine temperature of the sample 
from local observables as ksT = nS{0)/K.. 

5. Conclusion 

We demonstrated the extraction of density-density correlations and static structure 
factors from in situ images of 2D Bose gases. Careful analysis and modeling of the 
imaging response function allow us to fully eliminate the systematic effect of imaging 
imperfections on our measurements of density-density correlations. For thermal gases, 
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our measurement of the static structure factor agrees well with theory. For interacting 
2D gases below the BKT critical temperature, intriguingly, we observe positive density- 
density correlations in weakly interacting samples (g ^ 1) and negative correlations in 
the strong interaction regime (g'eff = 1-0). For all interacting gases, our static structure 
factor measurements agree with the prediction of the FDT as S{0) = ksTK/n. Extension 
of our 2D measurement can further test the prediction of anomalous density fluctuations 
in a condensate [23l HOj SU [42] and strong correlations in the quantum critical region 
j43l [44] . Finally, our analysis can be applied to perform precise local thermometry [45j 
and can potentially be used to extract the local excitation energy spectrum through the 
application of the generalized fluctuation-dissipation theorem j25i l28] . 
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Appendix A. Full analysis of the point spread function and the modulation 
transfer function 

Appendix A. 1. Point spread function 

Here, we describe our approach to characterize imaging imperfections using extended 
Nijboer-Zernike diflFraction theory [46j. To obtain the point spread function from the 
exit pupil function p{rp^9p)^ it is convenient to first decompose the exit pupil function 
using a complete set of orthogonal functions on the unit disk in the polar coordinates 

oo n 
n=Om=-n ^ 

where Z^{p^9) = B}^\p)e^^^ is a Zernike polynomial, the radial function R^{p) = 
Eto"^^' ki[in+Jj2-m^^^^ terminates at p = 1, and i?;^^ = when n - m is 

odd. The expansion coeflicient is given by 

'Ti _|_ 1 ra r27T 

/5r = ^/ / p{rp,ep)Z-^{rp/a,ep)rpdrpdep. (A.2) 
Tra^ Jo Jo 

If we then apply the expansion Eq. ( |A.lD to the Fourier transform of the exit pupil 
function p{k^9) = Jq Jq^ p{rp^9p)e~'^^^p^^^^^~^p'>rpdrpd9p and carry out the integration 
using the Zernike-Bessel relation Jq K^{p)Jjn{^p)pdp = (— we arrive 
at the following formula 

p{k,9) = 27ia^f2 E ^-/3-e^-(^+^)^^^^^M, (A.3) 

n=0 m=-n 
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Figure Al. Analysis of the imaging aberrations and the point spread function, (a) 
shows the wavefront aberration (the phase of the exit pupil function) determined from 
the fit to the experiment, (b) shows the expansion coefficients determined from 
Eq. (A.2), using the exit pupil function in (a), (c) shows a line-cut of the derived 
point spread function (solid line). The unaberrated point spread function is plotted 
for comparison (dashed line). Inset shows the 2D distribution of the aberrated point 
spread function and red line indicates the direction of the line-cut. Image size is 
(33 /im)2. 



where Jn{z) is the n-th order Bessel function of the first kind. The point spread function 
V{r^9) is the real part of e^^^p{k^9)\k=r/ad with proper normalization 

V{r,d) = ^Y. E (A.4) 

where M = 27Td'^^[e'^^^p{rp)]\rp=o = 27rrf^cos55 is the normalizing factor such that 
/ V{Y)d^r = 1. For a unaberrated system, only (3^ is non-zero and the above equation 
reduces to the form J\[z)jz^ as expected. 

Using the above equations, we can derive the point spread function from the fitted 
exit pupil function shown in Fig. |Al|(a). We calculate the expansion coefficients and 



evaluate the corresponding point spread function, see Fig. Al(b-c). 



Appendix A.2. Modulation transfer function 

It is straightforward to evaluate the modulation transfer function A^(k) = |P(k)| 
and the imaging response function A^^(k) = |P(k)p directly from the exit pupil 
function p(rp, Q^. Since the point spread function can be written as V{y) = [e'^^^p(k) + 
e~'^^^p*(k)]/4:7Ta'^Af\\^=r/adi its Fourier transform is 

V{k) = e + 7r) + e-V(r„ ^)]|.,.,«,, (A.5) 

where k = |k| is the spatial frequency and 9 is the polar angle of k in the image 
plane. From Eq. (A.5), the imaging response function is A^^(k) oc \p{kad^9 + tt) + 



e ^^'^^]9*(/carf, This result shows that the phase shift Sg is important since A^^(k) 
depends on the interference between p{kad^ O + n) and p^'^kad^ 9). The transmittance [/, 



defined in the exit pupil function Eq. (10), accounts for the radial envelope in (k). 



leading to the sharp edge at = ^ = 27rN.A./A. Either the continuous function 
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Eq. (A. 5) or the discrete Fourier model Eq. (12) can be used to calculate the imaging 
response function. 
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